Excitation of multiple wakefields by short laser pulses in dense plasmas 
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We present a theoretical investigation of the excitation of multiple electrostatic wakefields by the 
ponderomotive force of a short electromagnetic pulse propagating through a dense plasma. It is 
found that the inclusion of the quantum statistical pressure and quantum electron tunneling effects 
can qualitatively change the classical behavior of the wakefield. In addition to the well known plasma 
oscillation wakefield, with a wavelength of the order of the electron skin depth (Ae = c/cope, which 
in a dense plasma is of the order of several nanometers, where c is the speed of light in vacuum 
and Ljpe is the electron plasma frequency), wakefields in dense plasmas with a shorter wavelength 
(in comparison with Ae) are also excited. The wakefields can trap electrons and accelerate them to 
extremely high energies over nanoscales. 

PACS numbers: 52.35.Fp, 52.35.Mw, 52.38.Kd 
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Three decades ago, Tajima and Dawson [l| demon- 
strated that intense laser pulses can efficiently generate 
electron plasma waves (EPWs) in their wake as they 
travel through a low density plasma. Physically, the 
ponderomotive force of an intense laser pulse pushes 
electrons locally in plasmas, which in turn oscillate at 
the electron plasma frequency with respect to the neu- 
tralizing background of immobile positive ions. The dis- 
placement of the electrons within the plasma gives rise to 
electric fields, which can be much larger than any fields 
possible in a non-ionized material. An electron beam 
can surf on the electric field of a plasma wave picking 
up energy from the EPWs just as a surfer picks up en- 
ergy from a water wave in the ocean The idea of 
Tajima and Dawson has now been experimentally veri- 
fied SSiSii- 

Recently, there has been a great deal of interest in in- 
vestigating the properties of high-energy density plasmas 
that are created by high intensity laser pulses. To probe 
dense matter, such as those in the interior of white dwarf 
stars and Jovian planets, powerful laser-produced x ray 
sources have been developed. They produce monoener- 
getic line radiation capable of penetrating through dense 
and compressed materials at solid densities and above 
[lo| . The x-ray measurement techniques are indicative of 
a dense Fermi-degenerate plasma state in laboratories. In 
a dense Fermi plasma, the electron degeneracy leads to 
a consideration of the Fermi-Dirac electron distribution 
and electron quantum tunneling through the quantum 
Bohm potential lul, ll^, |l3| . Furthermore, there are also 
the spin force and spin magnetized electron current due 
to electron- 1/2 spin in dense magnetoplasmas 14, ITst . 



electron dynamics, and therefore one encounters numer- 
ous novel collective interactions in dense quantum plas- 
mas. Specifically, it should be stressed that the quantum 
Bohm force effect, arising from the finite width of the 
electron wave function, gives rise to dispersion of EPWs 
at nanoscales, which has important consequences to lo- 
calized EPW structures [l3l and plasmonic turbu- 
lence 
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The quantum statistical pressure, the quantum Bohm 
force and the quantum spin force drastically affect the 



In this Letter, we present a theoretical investigation 
of the multiple EPW (wakefield) excitation by the pon- 
deromotive force of a short laser pulse, accounting for the 
quantum statistical pressure and quantum Bohm force ef- 
fects in the EPW dynamics. For our purposes, we shall 
use the quantum Madelung fluid equations [I4I, which 
are composed of the electron continuity and electron mo- 
mentum equations, together with the Poisson equation, 
and derive the EPW (wakefield) equation in the presence 
of the radiation pressure. Chosing a specific form for the 
laser envelope, we solve the EPW equation analytically 
and numerically. We find that, due to the quantum Bohm 
force, multiple wakefields at nanoscales appear in dense 
quantum plasmas. In addition to the well known wake- 
field, with a wavelength of the order of the electron skin 
depth, a short wavelength wakefield is also excited, with 
a scale length approaching the Compton wavelength. It 
turns out that for the case of a laser pulse in the optical 
regime, the short scale wakefield is suppressed. However, 
for short laser pulse lengths and/or high density plas- 
mas, the energy density of the short scale wakefield may 
be comparable to that of the long wavelength wakefield. 
The consequences of our results are discussed. 

We consider the propagation of a high-frequency laser 
pulse, with the vector potential A = A exp{ikx — icut) 
-fee, in an unmagnetized dense plasma. Here c.c. stands 
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for the complex conjugate. The ponderomotive force of 
the high-frequency laser pulse drives longitudinal EPWs 
(wakefields) with a frequency much smaller than oj, but 
fast enough for the dynamics to take place on the elec- 
tron timescale. The ions form a neutralizing background 
in our dense plasma. The governing equations for the 
wakefields are then the electron continuity equation 



dnx 
~dt 



V • (nov) = 0, 
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the electron momentum equation 



at m 
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and the Poisson equation 
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rVV^m, (2) 
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A simple special case is found when Vg 
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results in exponentially damped wakefields 
below we will consider the case where we can have mul- 
tiple oscillatory wakefields, which occur for plasmas of 
moderate density (see the condition below). We start by 
investigating the solutions to the left-hand side of Eq. 
(5) in the absence of the driving laser field. Making the 
ansatz $ cc exp(ifc^), we obtain ^ — k^k'^ + kf = Q, where 
kl = - Vp) If? and kl = Am^u^^jt? . This has 

the solution 
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where ni is the electron density perturbation in the equi- 
librium value no, v is the electron fluid velocity per- 
turbation, $ is the Wakefield potential, e is the magni- 
tude of the electron charge, m is the electron rest mass, 
Vf — (27r/i/V3m,)(3no/87r)^/'^ is the Fermi speed, and h 
is the Planck constant divided by 27r. We have thus as- 
sumed that the Fermi electron pressure dominates over 
the electron thermal pressure, appropriate for a high den- 
sity plasma of moderate or low electron temperature. 
Several comments are in order. The second term in the 
right-hand side of (2), which represents the light pres- 
sure or the laser ponderomotive force, comes from aver- 
aging (over the laser period) the advection and nonlinear 
Lorentz force in the electron equation of motion Q . The 
fourth term in the right-hand side of (2) is the quantum 
Bohm force involving quantum electron tunneling in a 
dense quantum plasma. 

Combining Eqs. (l)-(3) we obtain the plasma wakefield 
equation in the presence of the light pressure in our dense 
plasma 
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where uipe = {A-KriQe^ /rrieY^'^ is the electron plasma fre- 
quency. In dense laboratory plasmas and in compact 
astrophysical objects, the latter is in the x-ray regime. 

Let us now consider the excitation of a one-dimensional 
wakefield by the high-frequency laser pulse that is prop- 
agating with the group velocity. Thus, we look for sta- 
tionary solutions of (4) in a comoving frame. Letting 
^ — X — Vgt, where Vg is the group velocity, we obtain 
from Eq. (g]) 



Unless the plasma density is very high (or w 



(3) 3> A:^, and the two solutions separate into oscil- 



lations with very different scale lengths. Denoting the 
solutions with k+ and fc_, we have for A:^/4 ^ 
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The analysis below will be valid for arbitrary values of 
and fc_, however, as long as both modes are oscilla- 
tory, which holds whenever Vg > vf and m |wg — > 
hujp apply. When the second inequality is a strong one, 
the oscillations at fc+ is close to the Compton scale (at 
least when Vg — Vp c^), and the fc_ oscillation is the 
standard plasma oscillation wakefield, although some- 
what modified by the inclusion of the Fermi pressure. As 
a starting point, however, no separation in magnitude of 
the two scales will be assumed. Firstly, we rewrite Eq. 
dSD as 



where 77 = 2muj^/f?c^. Eq 
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$ = ,7|Ap, (6) 
can be integrated to yield 
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where we have assumed as boundary condition that $ is 
zero before the arrival of the high-frequency laser pulse. 

Next, we assume that the laser pulse profile is a Gaus- 
sian, viz. |yl|2 = exp(— ^'2/L2)^ The energy density 
of the high-frequency laser pulse is assumed to be high 
enough such that the high-frequency laser pulse is chang- 
ing its shape on a longer timescale, as compared to the 
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Wakefield generation process. The simplest case is that 
of a very short pulse, i.e. when L <C kZ^, when the 
Wakefield after the laser pulse passage can be written as 
$ — <f>+ sin(fc+^ + (p+) + $_ sin(fc_^ + ip-), where (p+ 
and if- are constant phase angles and the amplitudes 
of the wakefields are proportional to the high-frequency 
laser pulse energy, and given by 

/oo 

In general, the wakefield amplitudes after the laser 
pulse passage can be written as 

such that the electric field amplitudes E± of the different 
wakefields in the case of a Gaussian profile is given by 

E± - (fc2 !fe2) W^exp(-4LV4) (8) 

where we have introduced W — \-^\^{£.')dC pro- 
portional to the high-frequency laser pulse energy. 
(Moreover, we have used cos{y) exp{—y^ / L^)dy = 
LV^e-^'/4) Thus, we see that efficient wakefield gen- 
eration requires that the laser pulse length is not much 
longer than the respective wakefield wavelength. If the 
Fermi speed is much smaller than the speed of light in 
vacuum, and lo ^ ujp, significant excitation of the long 
wavelength wakefield requires that the laser pulse length 
is not much longer than the electron skin depth, whereas 
significant excitation of the short scale wakefield requires 
a laser pulse length not much longer than the Compton 
wavelength. This condition of an extremely short pulse 
can be relaxed to some extent in a high-density plasma 
where the Fermi speed approaches the speed of light in 
vacuum, such as in a white dwarf star. For such densities 
we may have significant excitation also of the short scale 
wakefield for laser pulse lengths L ^ h/m{vg — VpY^^ or 
shorter. The full profile of the density perturbation ni 
induced by the high-frequency laser pulse is calculated 
numerically from Eqs ^ and ([7]) for the plasma number 
density no ~ lO^^cm"'^ (which is achievable in inertial fu- 
sion experiments) and a pulse length L of a few Compton 
wavelengths (i.e. k^L ~ 7). The result is shown in Fig. 
1. Due to the large separation in length-scales, the quan- 
tum scale wakefield is hidden to some extent in the large 
scale plot. However, the quantum scale oscillations are 
clearly shown in the zoom made in the sub-panel of Fig 1. 
For a higher plasma number density hq ~ 10'^*'cm~^ (as 



relevant for e.g. white dwarf stars) and a slightly shorter 
pulse length, the quantum scale oscillations can be di- 
rectly seen, as pictured in Fig 2 (pulse length /c+L ~ 4) 
and Fig 3. (pulse length k+L o±3). 

To summarize, we have presented a theoretical study 
of the EPW excitation by short laser pulses in a dense 
plasma. For this purpose, we have used the electron con- 
tinuity and electron momentum equation including the 
quantum statistical pressure, the quantum Bohm force, 
and the light pressure, together with the Poisson equa- 
tion to derive the driven EPW equation. The latter 
shows that the quantum force produces an additional 
EPW dispersion, besides that coming from the quantum 
statistical pressure. Our results reveal that the quantum 
Bohm force provides the possibility of a short wavelength 
(in comparison with the electron skin depth) wakefield, 
in addition to the standard plasma oscillation wakefield, 
with a wavelength of the order of the electron skin depth. 
For current experiments with the laser pulses in the opti- 
cal regime, the short scale wakefield is suppressed. How- 
ever, for very short laser pulses, of the order of a few 
Compton wavelengths, or for high density plasmas where 
the Fermi speed approaches the speed of light in vacuum, 
the energy density of the short scale wakefield is compara- 
ble to that of the long wavelength wakefield. The energy 
loss of the laser pulse is then correspondingly higher. The 
laser excited wakefields can trap electrons and accelerate 
them to high energies at nanoscales in dense plasmas, 
such as those in the next generation intense laser-solid 
density plasma experiments and in compact astrophysi- 
cal objects [3| (e.g. the interior of white dwarf stars and 
the planet Jupiter). Finally, the present investigation can 
be readily generalized by including an external magnetic 
field, in which case the propagation of both laser pulses 
and wakefields will be greatly affected. 
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FIG. 1: The total normalized plasma density wakefield 
for no — lO^^cm"^ and a Gaussian laser pulse profile 
A§ exp(— ^^/L^)with fc+L ~ 7. A zoom is made in the sub- 
panel to make the quantum scale oscillations clear. The nor- 
malizations used are ^ = '^p£,/c and n = eorjA^^n/qe. 
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FIG. 2: The total normalized plasma density wakefield 
for no ~ 10^°cm~^ and a Gaussian laser pulse profile 
Agexp(-CV-f'^)with k+L ~ 2.1. 
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FIG. 3: The total normalized plasma density wakefield 
for no — 10^°cm~^ and a Gaussian laser pulse profile 
A§exp(-^V-t'^)with k+L ~ 1.6. 



